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We report state-of-the-art atomistic simulations combined 
with high-fidelity conductance calculations to probe structure- 
conductance relationships in Au-benzenedithiolate(BDT)-Au 
junctions under elongation. Our results demonstrate that large 
increases in conductance are associated with the formation of 
monatomic chains (MACs) of Au atoms directly connected to 
BDT. An analysis of the electronic structure of the simulated 
junctions reveals that enhancement in the .s-like states in Au 
MACs causes the increases in conductance. Other structures also 
result in increased conductance but are too short-lived to be de- 
tected in experiment, while MACs remain stable for long simula- 
tion times. Examinations of thermally evolved junctions with and 
without MACs show negligible overlap between conductance his- 
tograms, indicating that the increase in conductance is related to 
this unique structural change and not thermal fluctuation. These 
results, which provide an excellent explanation for a recently ob- 
served anomalous experimental result [Bruot et al, Nature Nan- 
otech., 2012, 7, 35-40], should aid in the development of mechan- 
ically responsive molecular electronic devices. 

While benzene- 1 ,4-dithiolate (BDT) has been widely studied over 
the years for applications in molecular electronicsji"'^ recent dis- 
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coveries * of its tunable conductance properties have generated a 
renewed interest in the molecule. These discoveries include coun- 
terintutive conductance increases (exceeding an order of magnitude) 
during elongation of a Au-BDT-Au junction," 5 -! and a wide conduc- 
tance window spanning three orders of magnitude. 6 The tunability 
of BDT's conductance is enabled by the proximity of its highest oc- 
cupied molecular orbital (HOMO) to the Au electrode Fermi level 
(S.F )• Relatively small increases in the HOMO level result in reso- 
nant (or near resonant) tunneling, which, as demonstrated by Bruot 
et a/Scan be achieved through mechanical elongation of a Au-BDT- 
Au junction. However, elongation does not guarantee increases in 
co nduct ance, as evidenced by numerous previous Au-BDT-Au stud- 
ies 1 2 6 where conductance increases were not reported. The reason 
for this discrepancy is that the exact location of the HOMO level de- 
pends on the structural evolution of a junction, which may be influ- 
enced by the experimental setup and/or conditions. In order to make 
use of the desirable properties of BDT, and to facilitate experimental 
reproducibility, it is essential to determine the structure(s) responsi- 
ble for increasing conductance. 
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Prior computational work™ has thus far provided considerable 
insight into the behavior of BDT and the underlying mechanisms 
that control conductance. However, most prior computational stud- 
ies of molecular junctions have adopted simplifications which may 
ultimately make it difficult to connect the predicted behavior with 
experiment. For instance, in typical computational studies of BDT, 
the molecule is sandwiched between ideal, planar surfaces with an 
arbitrary initial geometry, and then stretched via geometry optimiza- 
tions. 
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is approach may not adequately capture many impor- 
tant aspects found in experiment, such as temperature effects, elon- 
gation rate effects, and non-ideal tip geometry. In order to better con- 
nect theory to experiment and facilitate a deeper understanding of 
the structure-conductance relationship, a different computational ap- 
proach is needed that can capture these environmental effects and the 
stochastic nature of junction formation. Here, we combine detailed 
atomistic simulations with first-principles calculations to investigate 
the structural basis for the conductance behavior of elongating Au- 
BDT-Au junctions. We perform a statistically significant number of 
simulations (>100) that, unlike previous comptuational studiesJ^Eni 
accounts for the spontaneous connection of a BDT molecule between 
two ruptured nanowire tips and thus are significantly less biased by 
their starting configuration and more representative of experiment. 5 



Specifically, in this work, we apply a hybrid molecular dynamics- 
Monte Carlo (MD-MC) protocolE^tMl within LAMMPS^and an in- 
house code to incorporate more realistic environmental features into 
computational studies of molecular junctions. Au-Au interactions are 
described using a semi-empirical many-body potential^ and the S- 
Au interactions are modeled with multi-site bonding curves parame- 
terized from density functional theory (DFT). — The simulation pro- 
cedure is illustrated in Fig. 1 and involves trapping a single BDT 
molecule between ruptured Au nanowire tips by mechanical deform- 
ing a BDT coated nanowire; once a BDT connects between two tips, 
the junction is evolved via mechanical elongation at finite tempera- 
ture (77 K), and the atomic coordinates are periodically extracted for 
use as input to electron transport calculations. This simulation proce- 
dure closely models the widel y used mechanically controllable break 
junction experimental methodEEHl and thus should produce experi- 
mentally representative configurations. Additionally, this procedure 
enables us to more efficiently collect a statistically relevant number 
of configurations than typically used in quantum mechanical simula- 
tions, allowing us to capture the ensemble behavior; specifically in 
this work, we perform 104 independent simulations to model break 
junction experiments, where we note that recent quantum mechan- 
ical studiesL^mi considered only at most four structurally distinct 
junctions. We believe that the geometries predicted by our simu- 
lations should better represent break junction experiments than pre- 
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Fig. 1 Simulation snapshots of the elongation of a BDT-coated Au 
nanowire, leading to the formation, elongation, and eventual rupture 
of a Au-BDT-Au junction. From top to bottom, Ax = 0.0, 4.0, 8.0, 
12.0, and 16.0 A. Monolayer molecules are removed after junction 
formation to isolate electrode geometry effects. All images in this 
work were rendered using Visual Molecular Dynamics.El 

vious computational workHSl due to the inclusion of (1) realistic 
electrodes, (2) finite temperature effects, (3) better statistics, and (4) 
junctions that are less biased by starting configuration. For more de- 
tails about the simulation method and applied force fields please refer 
to the ESI.t 

In our break junction simulations, a BDT molecule attaches be- 
tween two ruptured nanowire tips in 31 out of 104 (30%) indepen- 
dent runs, in excellent agreement with values reported in Au/BDT 
break junction experiments (30-40%). Akin to the experimental sit- 
uation, the formation of a molecular junction is dictated by the prox- 
imity of adsorbed BDT molecules to the nanowire fracture location, 
and the availability of bonding sites on the Au tip(s). The geometry 
of a Au tip depends on the structural pathway of the nanowire dur- 
ing the elongation process, and hence can vary substantially between 
runsJ 19 " 23 lThus, a large number of junction geometries are possible, 
which we expect to better represent experiment than prior computa- 
tional studies. 

To elucidate the structure-conductance relationship, we calcu- 
late the zero-bias conductance of our simulated junctions within the 
SIESTA 24 and SMEAGOlPSI packages, employing an approxi- 
mate self-interaction correction (ASIC)' 2 - 7 -'mefhod to the local density 
approximation (LDA), which yields conductance va lues that more 
closely match the experimentally most-probable value^I^SIfhan stan- 
dard DFT 29 (although they are still higher by a factor of approxi- 
mately six). The conductance values we report here for our simulated 
junctions are generally higher than the most-probable value reported 
in some experiments,' 2 ^ 2 ^' while they are within the range of values 
reported in more recent experiments.! 5 -!^ The source of the quantita- 
tive disagreement between the experimenally measured and theoreti- 



cally calculated values of Au-BDT-Au conductance is still a debated 
topic.' 30 " 32 ' As such, here we focus on relative changes in conduc- 
tance rather than the exact values. A detailed description of the con- 
ductance calculations, and results from benchmark calculations val- 
idating the conductance calculations and force fields are provided in 
the ESI.t 

In Fig. 2 we plot the conductance evolution of eight representative 
Au-BDT-Au junctions under elongation. Distinct differences in the 
shape of the conductance traces are apparent in Fig. 2; for clarity of 
presentation, we group similar traces into two separate plots shown 
in Fig. 2a and Fig. 2f. The curves in Fig. 2a undergo large gradual 
increases, while the curves in Fig. 2f are relatively flat; we note that 
both behaviors closely match recent experimental results 5 and that 
previous idealized computational studies have failed to capture this 
range of behaviors. 7 8 Furthermore, these results demonstrate that in- 
creases in conductance depend on structure and are not a natural con- 
sequence of the increasing potential energy of a junction (see Figure 
S5 in the ESIf). 

Interestingly, the breaking geometries shown in Fig. 2b-e all con- 
tain a monatomic chain (MAC) of Au atoms directly connected to 
BDT. Of the 3 1 simulations resulting in molecular junction forma- 
tion, 13% go on to form direct MAC -BDT connections during elon- 
gation. Here we find that each MAC-containing geometry in Fig. 2b- 
e results in conductance > 0.2Go- Closer examination of the entire 
stretching trajectories reveals that the initial jumps to conductance 
values larger than 0.1 Go coincide with the appearance of a MAC at 
a BDT-Au interface, suggesting that the BDT-MAC connection is re- 
sponsible for the increased conductance. To test this, we manually 
insert MACs at the BDT-Au interface of an ideal junction, as shown 
in Fig. 2k-m, finding that conductance increases as MACs are in- 
serted into otherwise static geometries. Analysis of the electronic 
structure of the simulated junctions reveals that MACs broadly (+/- 1 
eV) enhance the projected density of states (PDOS; see Figure S7 in 
the ESIf) around Zf for the MAC Au s and p z states, thus inducing a 
stronger coupling of the molecular states around Zf with the Au elec- 
trodes, resulting in a higher transmission. Importantly, this enhance- 
ment only occurs when BDT is bonded to a Au atom that bonds with 
one other Au atom; adatoms or atomically sharp tips (Fig. 2k) do not 
result in this enhancement (see Figure S7 in the ESIf). Previous stud- 
ies employing ideal, planar electrodes have reported increased con- 
ductance for BDT connected to a Au adatom at both electrodes.l^LU 
Here, we propose a distinct structural mechanism in which BDT is 
connected to a MAC at only one of the electrodes. This is supported 
by atomistic simulations that closely mimic break-junction experi- 
ments, providing evidence for the formation of MACs under realistic 
conditions. Meanwhile, prior studies employed manually generated 
and/or geometry-optimized (effectively K) structures, providing no 
evidence of the formation of such structures under dynamic, non- 
ideal conditions. 

In contrast to MACs, other low-coordination electrode structures 
such as Au-Au2-Au units (see Fig. 2i,j) do not increase conductance 
(note that Au-Au2-Au structu res ha ve been previously observed in 
DFT-based studies of Au tip s l 21 ' 33 !). Meanwhile, the structures in 
Fig. 2g,h result in high conductance, but their lifetimes are extremely 
short; in non-stretching MD simulations the structures shown in Fig. 
2g,h remain stable for less than 1 .0 ns, which is a time scale that is too 
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Fig. 2 Conductance behavior of Au-BDT-Au junctions undergoing elongation, (a) Four gradually increasing conductance traces, with the 
corresponding geometries immediately prior to junction rupture shown to the right. The individual traces are offset along the x-axis for clarity. 
At a distance of 0.1 A prior to junction formation the conductance is assumed to be IGo while at rupture the conductance is set to OGq. (f) 
Four relatively flat conductance traces, with the corresponding break geometries shown to the right, (k-m) The effect on conductance of 
manually connecting Au MACs (indicated with arrows) to BDT. The BDT is initially connected between two ideal Au(100) tips. A 
comparison of the transmission and projected density of states is included in the ESI, see Figures S4 and S5.t 



short to be measured in experiment. On the other hand, three of the 
four MACs in Fig. 2b-e remain stable for the complete duration of a 
1.0-/JS simulation without stretching (the actual simulated junctions 
were taken 0.2-0.5 A of elongation prior to the rupture structures 
shown in Fig. 2b-e), indicating that MACs may possess sufficient 
stability to be detected in experiments performed at 77 K. 

Note that MACs have been observed in experimentsii2M|35| an( j 
in simulation s ! 19 '' 21 ! f elongating Au nanowires without BDT. 
MACs have also been observed in simulations of thiolate-terminated 
molecules being pulled away from step edges on Au surfaces^! and 
in simulations of Au-alkanedithiolate-Au junctions. 37 Thus, the ap- 
pearance of MACs in Au/BDT break junction experiments is not 
without precedent and their formation is not an artifact of our simula- 
tion model. However, the significant impact of MACs on molecular 
conductance has not been previously demonstrated. 

To confirm that the conductance increases are not the result of the 
limited time scale of our simulations, we next compare the conduc- 
tance distributions of three structurally distinct Au-BDT-Au junc- 
tions (see Fig. 3) evolved at 77 K without stretching the junction. 
Each junction is evolved for 40 ns, with conductance measurements 
taken every 0.2 ns; this time frame allows for sufficient sampling 
of the most probable structures and results in converged histograms. 
The conductance of BDT connected to a MAC is statistically distinct 
from the two other deformed junctions with almost no overlap be- 
tween the histograms, indicating that the increases seen in Fig. 2a 
are indeed associated with their own unique structural motif and not 
simply a short-lived or improbable configuration. We also observe 
that the conductance fluctuations for BDT connected to a Au MAC 
is approximately double that of the other structures due to increased 
electrode motionM^J Increased fluctuations in BDT-MAC geometries 
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Fig. 3 Conductance histograms of three thermally evolving 
Au-BDT-Au junctions. The right-most junction/histogram 
corresponds to a system with a MAC, where we see negligible 
overlap with the non-MAC-containing junctions. The bin width is 
0.006Go- Standard deviations of the histograms (from left to right) 
are 0.014G , 0.018G , and 0.034G . 



have been reported previously, 38 albeit for unrealistic, manually ad- 
justed geometries, making it difficult to assess their connection to 



1-| | 3 



experiment. We also note that the histograms in Fig. 3 are similar to 
those in experiments by Tsutsui et al. ,1221 where a high-conductance 
state at 0.1 Go exhibited a larger peak width than the peak at lower 
values (0.01 G ). 

The similarities between our simulations and the experiments of 
Bruot et al .^suggest that MAC formation is the likely cause of the 
large (factor of five or more) gradual increases in conductance ob- 
served in these experiments. The lack of large conductance increases 
in experiments performed at 300 K 5 suggests that the structural mo- 
tif responsible for increased conductance has a distinct thermal de- 
pendence {i.e., the structure is unstable at higher temperature). It 
has been previously 20 established that the formation of MACs in 
mechanically deformed Au nanowires will depend on temperature, 
where MACs form with the highest frequency and stability at low 
temperature. We confirm this by performing 298-K simulations with- 
out stretching of the MAC geometries in Fig. 2b-e (the actual sim- 
ulated junctions were taken 0.2-0.5 A of elongation prior to the rup- 
ture structures shown in Fig. 2b-e). In all four cases the junction 
undergoes thermoactivated spontaneous breakdown in less than 1.0 
ns, indicating that MACs are not stable at high temperature. Thus, 
the thermal instability of MACs explains why large conductance in- 
creases were not observed in 300-K experiments. 5 While in the work 
of Bruot et al™ there were no large conductance increases at 300 
K, large conductance values were reported in the 300-K conductance 
histogram, which may have been caused by other significant confor- 
mational changes to the junction at 300 K. 

Other fa ctors, such as Au adatoms and strained S-Au bonds, have 
been shown lW I 11 ! 40 ! to increase conductance in idealized junctions, 
and may also play a role under certain conditions. For instance, 
adatom formation may be important in STM-based break junction 
experiments 2 where a planar surface serves as one of the electrodes. 
Meanwhile, the small (~0.001Go) reversible changes in conductance 
(with respect to junction compression and elongation) in the work of 
Bruot et a/.l^are likely caused by reversible structural change such 
as a st raine d S-Au bond. On-hollow bonding geometry,^l22l high tilt 
angles,^E2 and S atoms embedded in Au contacts^ have all been 
proposed to explain large values of conductance, but these explana- 
tions do not seem likely given the very large, gradual conductance 
increases observed in the recent work of Bruot et al™ or in our com- 
putational studies. Additionally, we note that since the elongated 
junctions are under tension, it makes it more probable for BDT to 
adopt an upright geometry with low tilt angle, bonded with a single 
Au atom at each electrode. 

It is also important to note that the formation of MACs in prior 
theoretical work focusing on flat platesHSlmay be suppressed due to 
the electrode geometry choice (i.e., flat surfaces). Flat surfaces could 
discriminate against the formation of MACs since the dynamic struc- 
tural fluxionality (i.e., lengthening and weakening of Au-Au bonds) 
of flat Au surfaces is lower than that of nanostructured Au surfaces.^ 
The use of DFT-based geometry optimizations,!^ where elongation 
rate and temperature effects are neglected, may also explain the ab- 
sence of MACs in prior work, since the appearance of MACs has 
been linked to both the elongation rate and temperature, explore 
this we perform a Au-BDT-Au elongation simulation using ideal, pla- 
nar Au(100) surfaces to model the electrodes (see Fig. 4). In this 
simulation a Au adatom is first pulled from the surface of one of 




Fig. 4 Comparison of the conductance (top) and structural (bottom) 
evolution of Au-BDT-Au junctions where the electrodes are 
modeled as (left) nanowire tips and (right) planar surfaces. 



the electrodes, which then evolves into a 2-atom MAC before the 
junction ruptures. This result supports the important role of temper- 
ature in MAC formation, as simulations employing geometry opti- 
mizationsl^Kilhave not led to MAC formation. In support of our pro- 
posed structural mechanism, adatom and MAC formation in our sim- 
ulations is limited to a single electrode, while prior workEl reported 
significant conductance increases only when BDT connected to an 
adatom at both electrodes. Although the flat-surface and nanowire 
junctions exhibit roughly the same ductiliy (i.e., elongation prior to 
rupture), there are significant differences between the two curves in 
Fig. 4. Both tend to increase during elongation, however the increase 
is much larger for the nanowire tips. For the planar electrodes the 
conductance increases by a factor of almost 3, while for the nanowire 
tips the conductance increases by over an order of magnitude. The 
changes in conductance are also more gradual for the planar elec- 
trodes, whereas the nanotip data exhibit sharp fluctuations owing to 
the increased electrode motionM^I 

In summary, by combining realistic simulations of molecular junc- 
tion formation and elongation with high-fidelity conductance calcula- 
tions, we have provided important new insight into the conductance 
behavior of Au-BDT-Au junctions. Namely, we showed that BDT 
connected directly to a MAC results in enhanced conductance, and 
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is caused by enhancements in the .? and p z density of states around 
Zf in Au MAC atoms. This result offers an excellent explanation for 
the large, anomalous conductance increases observed in Au-BDT-Au 
break junction experiments^ and may additionally explain the large 
transmission observed in experiments by Kim and co-workers.l^ To 
further support our conclusion, we show that BDT-MAC structures 
are stable for long simulation times performed at 77 K, and exhibit 
very little overlap in their conductance distributions with those of 
other realistic junction geometries, and thus should be detectable in 
experiment. These results may aid experimental reproducibility by 
prompting the development of new strategies for tailoring the struc- 
tures and properties of molecular junctions, in particular mechan- 
ically reponsive devices.EI Our findings also stress the importance 
of incorporating realistic environmental conditions into simulation 
models of molecular junctions for better resolving experimental un- 
certainty. 
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